Explanation of the discrepancy between the measured and atomistically calculated 

yield stresses in body-centered cubic metals 
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We propose a mesoscopic model that explains the factor of two to three discrepancy between 
experimentally measured yield stresses of BCC metals at low temperatures and typical Peierls 
stresses determined by atomistic simulations of isolated screw dislocations. The model involves a 
Frank-Read type source emitting dislocations that become pure screws at a certain distance from 
the source and, owing to their high Peierls stress, control its operation. However, due to the 
mutual interaction between emitted dislocations the group consisting of both non-screw and screw 
dislocations can move at an applied stress that is about a factor of two to three lower than the stress 
needed for the glide of individual screw dislocations. 
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INTRODUCTION 



It has been firmly established by many experimen- 
tal and theoretical studies performed in the last forty 
years that the plastic behavior of body-centered-cubic 
(BCC) metals is controlled by 1/2(111) screw disloca- 
tions the cores of which are non-planar (for reviews see 
0,1111 Si, 110) However, the only direct experimen- 
tal observation that suggests such core spreading is the 
high-resolution transmission electron microscopic (TEM) 
study of Sigle while the primary source of our un- 
derstanding of the dislocation core structure and related 
atomic-level aspects of the glide of 1/2(111) screw dislo- 
cations is computer simulation. Such calculations have 
been made using a broad variety of descriptions of in- 
teratomic forces, ranging from pair-potentials [fll Ii(tI Iii1| 
to density functional theory (DFT) based calculations 
0, 0, Q| and studies exnploying other quantum me- 
chanics based methods jlj, [IgI, |13 ■ 

The vast majority of atomistic studies of the core struc- 
ture and glide of 1/2(111) screw dislocations in BCC met- 
als were carried out using molecular statics techniques 
and thus they correspond to K. A problem encoun- 
tered universally in all the calculations of the critical 
resolved shear stress (CRSS), i.e. the Peierls stress, at 
which the screw dislocation starts to glide, is that it is 
by a factor of two to three larger than the CRSS ob- 
tained by extrapolating experimental measurements of 
the yield and flow stresses to K. The following are a 
few examples. Basinski et al. 0] measured the flow 
stress of potassium in the temperature range 1.5 K to 
30 K and extrapolated to K to get 0.002^ to 0.003/^ 
where /z = (Cu - C12 + C44)/3 is the (111){110} shear 
modulus and Cu, C12, C44 are elastic constants. Similar 
values were found by Pichl and Krystian jl9j |. The val- 
ues of the CRSS when the maximum resolved shear stress 
plane (MRSSP) is a {110} plane, calculated using a pair 
potential derived on the basis of the theory of weak pseu- 
dopotentials is 0.007^ to 0.009/x J^J. More recently, 
Woodward and Rao [ll| calculated the CRSS in molyb- 
denum using the many-body potentials derived from the 



generalized pseudopotential theory |21| and a DFT based 
method. When the MRSSP is a {110} plane, they found 
the CRSS to be between 0.018/i and 0.020/i. A simi- 
lar value of the CRSS, 0.019^, was found in calculations 
employing the tight-binding based bond-order potential 
for molybdenum l22l |. Experimental measurements 
of HoUang et al. [23, extrapolated to K, give for the 
CRSS in molybdenum 0.006/x. A similar problem was 
encountered by Wen and Ngan j^l] who used the Em- 
bedded Atom Method (EAM) potential for iron and the 
Nudged Elastic Band method to analyze the activation 
enthalpies for kink-pair nucleation on screw dislocations. 
The calculated yield stress at K was about 0.01 3u while 
the experimental values, reported by Aono et al. |2^ are 
0.005/i to 0.006/x. This ubiquitous higher value of the 
calculated CRSS, found independently of atomic interac- 
tions, suggests that the origin of this discrepancy cannot 
be sought on the atomic scale of the motion of individual 
dislocations but rather on mesoscopic scale where a large 
number of elastically interacting dislocations glide at the 
same time. In this context it should be noticed that the 
only atomistic simulation that predicts yield stress close 
to that measured experimentally considered a planar dis- 
location network of 1/2[111] and 1/2[111] screw disloca- 
tions with [001] screw junctions j2y|. Such a network 
moved in the (110) plane at the stress about 50% lower 
than the Peierls stress of an isolated screw dislocation. 

In-situ TEM observations of dislocation sources in 
BCC transition metals showed that in thin foils straight 
screw dislocations formed near the source and moved very 
slowly as a group (S, El IH HI HI HI HI . Hence they 
fully control the rate at which the source produces dis- 
locations. In the foils used in TEM the applied stresses 
are very low but a similar control of the sources by sessile 
screw dislocations can be expected in the bulk at stresses 
leading to the macroscopic yielding. However, at higher 
stresses dislocations move faster and do not become pure 
screws immediately after leaving the source but at a dis- 
tance from the source. Indeed, even in situ observations 
at higher stresses do not show straight screw dislocations 
emanating directly from the sources |34| . 
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In this paper, we propose a mesoscopic model involv- 
ing a Frank-Read type source [s^ emitting dislocations 
of generally mixed character that become pure screw dis- 
locations at a distance from the source and, owing to 
their high Peierls stress, control its operation. However, 
there are a number of non-screw dislocations between the 
screws and the source, which can move easily. These dis- 
locations exert a stress on the screw dislocations and this 
stress, together with the applied stress, act on the screw 
dislocations by the force equal to that needed to overcome 
the Peierls stress. Screw dislocations can then move at 
an applied CRSS that is about a factor of two to three 
lower than the CRSS needed for the glide of individual 
screw dislocations. 



MODEL OF DISLOCATION NUCLEATION AND 
MOTION 

Let us consider a Frank- Read source (see e.g. jl^) 
that produces dislocations in a BCC metal. It emits, as 
always, dislocation loops that have a mixed character and 
expand easily away from the source since their Peierls 
stress is low. However, at a certain distance from the 
source, a significant part of the expanding loop attains 
the screw orientation and becomes much more difficult to 
move owing to the very high Peierls stress of pure screws. 
The rest of the loop, having a mixed character, continues 
to expand which leads to further extension of the screw 
segments. As a result, the source becomes surrounded by 
arrays of slowly moving screw dislocations, as depicted 
schematically in Fig.^ Further operation of the source 
is hindered by their back stress and effectively controlled 
by the ability of the screw dislocations to glide. 
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FIG. 1: Schematic operation of a dislocation source in bee 
metals. The curved non-screw segments migrate away, leaving 
behind a new pair of screw dislocations. 

The operation of the source is driven by the applied 
stress, a a, which acts by the Peach-Koehler force 
a ah (per unit length) on the dislocation that bows out. 
This dislocation obviously has a mixed character. Let 
us consider now that there are Ns screw dislocations at 
distances Xi from the source and Nm dislocations, gen- 
erally of mixed character, positioned between the source 



and the screw dislocations. We approximate the latter 
as straight lines of the same orientation as the screws, 
positioned at distances from the source, but with a 
negligible Peierls stress compared to that of the screws. 
In the framework of the isotropic elastic theory of dislo- 
cations the condition for the source to operate is then 
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where t is the line tension of the emitted dislocations, 
b their Burgers vector, R the half-length of the source, 
fi the shear modulus, and a, (3 constants of the order of 
unity. The first term is the force arising from the line 
tension that pulls the dislocation back and the second 
and third terms are forces produced by the stress fields 
of screw and non-screw dislocations, respectively, present 
ahead and/or behind the source. In the following we 
neglect the interaction between dislocations ahead and 
behind the source as they are far apart. Moreover, the 
dislocation sources are frequently single-ended (see e.g. 
[s^jl^ Hence we analyze only dislocations ahead of the 
source, i.e. those towards which the source bows out. 

It should be noted here that the screw dislocations in 
the array ahead of the source are not pressed against any 
obstacle and thus they do not form a pile-up. Within the 
approximations defined above, the ith screw dislocation 
will move provided 



ub-^ 1 lib 

CTa + T^- > + 7^ 



I 



27r ^ — ' Xi 



2'Ka ^-^ Xi 
fc=i * 



lib 1 
yk 27ra Xi 



> (jp 



(2) 

where txp is the Peierls stress of screw dislocations. The 
second and third terms are stresses arising from screw 
and non-screw dislocations, respectively, and the fourth 
term is the stress arising from the dislocation associated 
with the source that is also treated as a straight line of 
the same type as all the other mixed dislocations. Since 
the Peierls stress of non-screw dislocations is negligible, 
the ^th non-screw dislocation can move provided 
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The meanings of individual terms are analogous to those 
in equation (j^J. 

Now, the question asked is how large stress, tJa, needs 
to be applied so that the screw dislocations can move so 
far away from the source that they either reach a sur- 
face or encounter dislocations of opposite sign from an- 
other source and annihilate. In both cases the source 
then keeps producing new dislocations indefinitely. In 
the former case these dislocations keep vanishing at the 
surface and the latter case leads to the propagation of 
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slip through the sample. In order to investigate the 
problem formulated above, we performed the following 
self-consistent simulations for certain fixed values of the 
Peierls stress, dp, and applied stress Ua- First, we choose 
a half-length of the source, R, and a distance from the 
source, ymax, beyond which the expanding loop always 
attains the screw character. The first mixed dislocation 
emitted by the source becomes screw when reaching the 
distance ymax and then moves to a distance xi, deter- 
mined by equation (j^J. Provided that the source can 
operate, i.e. the inequality (Q) is satisfied, another dis- 
location is emitted from the source. The position of this 
dislocation is determined by equation if it does not 
reach ymax and by equation ^ if it does. Subsequently, 
the position of the first dislocation, xi, is updated to sat- 
isfy equation Q , which allows also the second dislocation 
to move. In this way a new position of the first disloca- 
tion, xi, and the position of the second dislocation, ei- 
ther yi if smaller than ymax or X2 if larger than ymax, 
are found self-consistently. This self-consistent process 
is then repeated for the third, fourth, etc., dislocations 
until the source cannot emit a new dislocation, i.e. when 
inequality is no longer satisfied. The result of this 
calculation is the number of screw dislocations, Ng, and 
mixed dislocations, Nm, as well as their positions ahead 
of the source, when the source becomes blocked by the 
back-stress from all the emitted dislocations. The first 
screw dislocation is then at a position xi = Xmax and fur- 
ther operation of the source can proceed only if this screw 
dislocation is removed, as argued above. The source can 
then continue operating in a steady-state manner, pro- 
ducing a large number of dislocations that mediate the 
macroscopic plastic flow. 

RESULTS 

In the following numerical simulations the applied 
stress, (Ta, has been set equal to O.Sctp and Q.bap, re- 
spectively, in order to investigate whether the source can 
operate at stress levels corresponding to experimental 
yield stresses extrapolated to K, as discussed in the 
Introduction. Three values of the Peierls stress, cp, have 
been considered that fall into the range found in atom- 
istic studies of transition metals [ij, [2^ |2J| , namely 
0.01/i, 0.02/i, and 0.03/x. Three different positions at 
which mixed dislocations transform into screw ones have 
been considered, namely ymax/b = 500, 1000, and 2000. 
The dependence on the size of the source, i?, was also in- 
vestigated. However, this dependence is very weak since 
R enters only through the line tension term in ^ and 
this is always small compared to the terms arising from 
the back-stress of emitted dislocations. Hence, without 
the loss of generality, we set R — ymax- The values of 
parameters a and /?, entering equations to Q have 
all been set to one and the usual approximation for the 



line tension, r = /i6^/2 [s^, adopted. 
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FIG. 2: Positions of mixed (open circles) and screw (full cir- 
cles) dislocations and the corresponding stresses when the 
source is blocked by the back-stress for the case Oajop — 0.5, 
o-p/^ = 0.02, and ymax/b = 500. 

Results of such simulation are presented in detail for 
(Ta = 0.5crp, (Tp = 0.02/i and ymax — 5006 in Fig.|2 
where positions of the dislocations ahead of the source 
and stresses acting on them are shown. In this case 
Xmax/ymax = 1-8. It should be uotcd that the stress 
exerted on the majority of screw dislocations is practi- 
cally equal to their Peierls stress. The distances Xmax 
found for the above-mentioned two values of CTq, three 
values of dp and three values of ymax, are summarized 
in TablcU 
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TABLE I: The distance which the leading screw dislocation 
advances from the source, Xmax, as a function of the distance 
ymax from the source at which dislocations become screw, for 
the applied stress Oajop and the Peierls stress of the screw 
dislocations crp/^. 

These results suggest that, for a given applied stress, 
the ratio Xmax/ymax is almost constant, independent of 
ymax, and only weakly dependent on the magnitude of 
the Peierls stress crp. At (la/crp = 0.3, most of the dislo- 
cations are mixed and Xmax/ymax ~ 1-3. With increas- 
ing stress, more emitted dislocations become screw and, 
at aa/crp = 0.5, Xmax/ymax ~ 2, which implies that 
the numbers of mixed and screw dislocations ahead of 
the source are very similar. Very importantly, the stress 
exerted on most of the screw dislocations is practically 
equal to their Peierls stress, see Fig.|21 
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CONCLUSION 



The distinguishing characteristic of the model pre- 
sented in this paper is that it does not consider the ghde 
of a single screw dislocation but movement of a large 
group of dislocations produced by a Frank-Read type 
source. In general, this source produces dislocation loops 
of mixed character that transform into pure screws at a 
distance ymax from the source. Hence, the group of dislo- 
cations consists of screw dislocations at distances larger 
than i/max and non-screw dislocations near the source. 
It is then the combination of the applied stress with the 
stress produced by the dislocations in this group that acts 
on the screw dislocations and is practically equal to their 
Peierls stress. However, after emitting a certain number 
of dislocations the source becomes blocked by their back- 
stress and, at this point, the leading screw dislocation 
reaches the distance Xmax from the source. Nonetheless, 
the source can continue operating if a dislocation of oppo- 
site sign, originating from another source, annihilates the 
leading screw dislocation. This requires an average sepa- 
ration of sources about 2xmax- Since the pinning points 
of the sources for a given slip system are produced by in- 
tersections with dislocations in other slip systems, their 
separation is related to the dislocation density in these 
systems. For example, in a deformed molybdenum crys- 
tal this density is of the order of 10^^ m~^ which im- 
plies separation of dislocations between 30006 and 40006, 
for the lattice parameter of Mo equal to 3.15 A. These 
values are in the range of 2xmax for applied stresses that 
are between one-third and one-half of the atomistically 
calculated Peierls stress for the sources of the size com- 
patible with the above-mentioned density of dislocations. 

The implication of the present study is that the val- 
ues of the Peierls stress of screw dislocations in BCC 
metals found in atomistic studies cannot be compared 
directly with the yield stress obtained by extrapolating 
experimental measurements to K. The experiments do 
not determine the stress needed for the glide of individ- 
ual screw dislocations but, instead, the stress needed for 
the operation of sources that are hindered by the sessile 
screw dislocations. These sources can operate at stresses 
lower than the Peierls stress owing to the collective mo- 
tion of screw and mixed dislocations produced by them, 
as described in this paper. Consequently, the discrepancy 
between the calculated Peierls stress and the measured 
yield stress is not a consequence of the inadequacy of 
the description of atomic interactions, which has often 
been raised as a possible explanation, but incorrectness 
of their direct comparison. 
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